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LONG-TERM  GOAL 

Our  long-term  goal  is  to  improve  the  aeeuraey  of  numerieal  predietion  models  of  weather  and  elimate. 
We  shall  eoneentrate  on  sourees  of  predietion  error  involving  interaetions  between  physieal  phenomena 
having  different  timeseales. 

OBJECTIVES 

We  examine  how  inadequate  representation  of  rapidly-varying  (i.e.,  stoehastie)  atmospherie  effeets  in 
General  Cireulation  Models  (GCMs)  ean  systematieally  affeet  predietion  on  a  variety  of  seales,  and 
how  errors  arising  from  this  inadequaey  may  be  ameliorated.  Coneurrently,  we  develop  useful,  effieient 
and  aeeurate  methods  of  aeeounting  for  these  systematie  effeets  of  stoehastie  foreing,  and  also  provide 
methods  for  estimating  the  spread  of  ensemble  predietions  taking  into  aeeount  multiplieative  stoehastie 
effeets. 

APPROACH 

We  approaeh  the  problem  both  analytieally  and  numerieally.  Linear  empirieal  models  and  barotropie 
models  of  the  atmosphere  present  eases  that  are  meteorologieally  relevant,  yet  simple  enough  to  solve 
analytieally,  and  we  use  these  models  to  identify  and  solve  some  of  the  numerieal  problems  that  will 
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arise  in  realistic,  more  complex  models.  We  intend  to  apply  our  results  to  studying  the  impact  of  sto¬ 
chastic  forcing  in  more  complete  models,  culminating  with  long  runs  of  stochastically  forced,  full- 
fledged  climate  GCMs  and  numerical  weather  prediction  (NWP)  models. 

The  theory  of  stochastic  white  noise  is  well  known  and  analytic  solutions  are  available,  at  least  in  prin¬ 
ciple.  We  have  extended  these  results  to  the  more  realistic  problem  of  stochastic  colored  noise  forcing, 
and  are  currently  investigating  how  best  to  implement  these  results  in  a  real-time  forecasting  system. 

Investigations  into  the  numerical  implementation  of  stochastic  theory  is  still  a  major  effort.  We  are  in 
communication  with  numerical  analysts  who,  at  our  request,  consider  the  errors  incurred  by  including 
stochastic  forcing  into  existing  deterministic  models,  most  of  which  use  a  leap-frog  integration  scheme. 
Further,  we  consider  recent  developments  in  solving  numerically  the  Fokker-Planck  equation  for  the 
forecast  distribution. 

WORK  COMPLETED 

We  have  completed  numerical  work  on  the  steadily  forced  barotropic  vorticity  equation  with  stochastic 
variations  in  the  basic  state.  We  also  completed  analysis  of  a  similar  system  with  stochastically  varying 
friction  on  the  forced  waves. 

We  derived  analytic  expressions  for  the  expected  mean  response  to  colored-noise  forcing  with  widely 
different  timescales. 

We  produced  linear  inverse  models  for  both  winter  and  summer.  These  models  produce  forecasts  of 
weekly  averaged  global  atmospheric  circulation  and  diabatic  heating  anomalies  whose  skill  is  competi¬ 
tive  with  NCEP’s  MRF  ensemble  at  week  2,  and  possibly  better  at  weeks  3  and  4.  We  have  analyzed  the 
models  and  determined  predictability  limits  of  the  models. 

RESULTS 

We  have  obtained  two  principal  results  this  year: 

The  fact  that  stochastic  variations  often  do  not  decay  on  timescales  much  faster  than  those  of  the 
weather  systems  to  be  predicted  required  us  to  consider  how  robust  are  the  white  noise  approximations 
used  in  our  previous  year’s  research.  We  developed  a  general  formula  for  the  drift  in  the  mean  response 
of  a  steadily  forced  linear  system  due  to  stochastic  perturbations  with  timescales  ranging  from  much 
faster  than  the  deterministic  timescales  to  much  slower  than  the  deterministic  timescales.  We  showed 
that  this  general  formula  is  rigorously  correct,  and  numerical  simulations  verified  this.  This  was  an 
important  advance  upon  analytic  results  last  year,  which  only  exactly  agreed  with  numerical  simulations 
at  asymptotically  long  and  short  decorrelation  time  scales  of  the  stochastic  variations.  We  found  that 
whether  stochastic  variations  in  the  basic  state  of  a  steadily  forced  barotropic  vorticity  model  act  to  sta¬ 
bilize  or  destabilize  Rossby  waves  depends  subtly  upon  their  impact  on  the  wave  properties  of  the  sys¬ 
tem. 

Secondly,  we  constructed  a  linear  inverse  model  (LIM)  suitable  for  studies  of  atmospheric  extratropical 
predictability  on  weekly  time  scales  using  global  observations  of  the  past  30  years.  Notably,  the  LIM 
includes  tropical  diabatic  heating  as  an  evolving  model  variable  rather  than  as  an  externally  specified 


forcing,  and  also  includes,  in  effect,  the  feedback  of  the  extratropical  weather  systems  on  the  more 
slowly  varying  circulation.  We  have  shown  that  both  of  these  features  are  important  contributors  to  the 
realism  of  this  empirical-dynamical  model. 


Predictability  of  Weekly  Averages  during  Winter 


(c)  Potential  Predictability 
with  Tropical  Heating 


7  days  14  days  21  days 


(d)  Potential  Predictability 
without  Tropical  Heating 


28  days  35  days  42  days 


1.  Forecast  skill  and  predictability  of  weekly  averages  during  winter. 

Top:  Correlation  of  observed  and  3-week  forecasts  of  upper  troposphere  streamfunction 
anomalies,  averages  over  52  forecast  cases  in  the  winters  of  1985/86 — 1988/89.  (a)  LIM. 

(b)  GCM  (N CEP  Reanalysis  model). 

Bottom:  Potential  predictability  limit:  forecast  lead  at  which  skill  (i.e.,  the  correlation  of 
observed  and  predicted  anomalies)  drops  below  0.5.  (a)  Determined  from  the  full  LIM. 

(b)  Determined  from  a  version  of  the  LIM  in  which  the  effects  of  tropical  forcing  are 
removed. 

At  week  2  the  LIM’s  foreeast  skill  is  eompetitive  with  that  of  the  global  nonlinear  medium-range  fore¬ 
east  (MRF)  model  with  O(IO^)  degrees  of  freedom  in  use  at  the  National  Centers  for  Environmental 
Predietion  (NCEP).  Figure  1  (a  and  b)  shows  sueh  a  eomparison  for  Week  3  foreeasts  made  throughout 


the  winters  of  1985/86-1988/89.  Analysis  of  the  LIM  suggests  that  without  inelusion  of  tropieal  heat¬ 
ing,  weekly  averages  in  the  extratropies  may  be  predietable  only  about  two  weeks  ahead,  but  with  trop¬ 
ieal  heating  ineluded,  they  may  be  predietable  as  far  as  seven  weeks  ahead.  This  diflferenee  is  displayed 
in  Figs,  le  and  Id.  This  shows  that  aeeurate  predietion  of  tropieal  diabatie  heating  is  key  to  enhaneing 
extratropieal  predietability 

IMPACT/APPLICATION 

Our  researeh  has  a  direet  impaet  on  ensemble  weather  and  elimate  predietion.  We  have  found  that  we 
needed  to  use  numerieal  methods  and  eriteria  for  ehoiee  of  time  step  and  integration  time  that  are  gener¬ 
ally  not  used  in  the  numerieal  modeling  eommunity.  The  methods  eurrently  used  in  stoehastie  numeri¬ 
eal  models  are  likely  to  be  highly  inaeeurate,  and  our  work  helps  to  1)  quantify  how  inaeeurate  those 
methods  may  be  and  2)  offer  alternatives  to  them  in  order  to  improve  the  aeeuraey  of  numerieal  predie¬ 
tion  models. 

Our  LIM  results  mark  the  first  time  that  a  sueeessful  linear  model  of  atmospherie  variability  on  weekly 
time  seales  has  been  eonstrueted.  If  sueh  variability  is  effectively  linear  then  a  detailed  knowledge  of 
nonlinearities  may  not  result  in  signifieant  further  improvement  in  predietability  on  these  time  seales. 

The  LIM  should  also  prove  useful  in  our  planned  researeh  this  year.  Beeause  it  produees  realistie  atmo¬ 
spherie  variability  yet  has  only  37  degrees  of  freedom,  it  is  partieularly  suited  to  analyzing  the  effeets  of 
realistie  multiplieative  noise.  We  have  already  shown  (see  last  year’s  report)  that  aoeurately  determining 
the  effeets  of  multiplieative  noise  takes  very  lengthy  eomputer  runs,  and  that  additionally  there  may  be 
unexpeeted  numerieal  modeling  pitfalls  in  the  eorreet  treatment  of  the  noise.  Although  we  have  learned 
mueh  about  the  effeets  of  multiplieative  noise  in  the  barotropie  eontext,  it  may  be  quite  a  leap  between 
that  and  a  full  GCM,  where  very  long  integrations  eould  be  impraetieal  without  additional  guidanee. 
Thus,  for  our  multiplieative  noise  experiments  the  LIM  may  be  a  very  important  intermediate  model 
between  the  barotropie  model  we  used  next  year  and  a  full  GCM  we  will  use  next  year. 

TRANSITIONS 

The  analytieal  work  we  have  done  has  proven  to  be  interesting  to  the  general  stoehastie  modeling  eom¬ 
munity.  We  diseussed  our  work  at  a  workshop  hosted  by  the  Centre  Europeen  de  Caleul  Atomique  et 
Moleeulaire,  in  Lyon,  Franee,  last  summer  and  have  been  invited  to  partieipate  in  a  workshop  on  sto- 
ehastie  elimate  modeling  hosted  next  June  by  the  Bereieh  Stoehastik,  Institut  fuer  Mathematik  der 
Humboldt  Universitaet  in  Chorin,  Germany.  Dr.  Sardeshmukh  has  presented  these  results  in  a  series  of 
invited  talks  at  the  EGS  meeting  in  Niee,  the  Geophysieal  Fluid  Dynamies  Laboratory  (GFDL)  in  Prin- 
eeton,  and  NCAR  in  Boulder. 

We  have  also  interested  the  numerieal  analysis  eommunity  and  Dr.  Penland  will  give  a  talk  in  Oetober  at 
the  Dept,  of  Applied  Mathematies,  University  of  Indiana,  Bloomington. 


We  have  used  the  LIM  to  produee  realtime  foreeasts  during  both  winter  and  summer.  These  foreeasts 
are  made  available  worldwide  at  the  CDC  web  site  (http://www.ede.noaa.gov/~erw/lim). 
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The  methods  developed  in  this  projeet  are  also  being  applied  to  examine  feedbaeks  between  the  oeean 
and  atmosphere  in  the  extratropies  on  submonthly  time  seales.  Dr.  Newman  is  eollaborating  with 
Miehael  Alexander  (CDC)  and  Gabriel  Lau  (GFDL)  on  this  projeet. 

Dr.  Sardeshmukh  and  Dr.  Penland  are  eollaborating  with  Gilbert  Compo  (CDC)  on  a  NOAA-funded 
projeet  involving  the  National  Centers  for  Environmental  Predietion’s  Medium-Range  Foreeast  Model. 
Using  a  distributional  approaeh,  they  evaluate  that  model’s  ability  to  prediet  extreme  temperature  and 
preeipitation  events  during  El  Nino,  La  Nina  and  normal  epoehs. 
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